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Abstract 

We present a method to find the number of real and imaginary observ- 
able parameters coming from the Yukawa sector in an arbitrary gauge theory. 
The method leads naturally to a classification of Yukawa couplings accord- 
ing to their symmetries and suggests a new parametrization of masses and 
mixings that is useful to study the behaviour of Yukawa couplings under the 
renormalization group. We apply it to some examples based on the Standard 
Model with Yukawa couplings obeying various chiral symmetries. We also 
show how our method of parameter counting can be used in some models with 
an enlarged leptonic sector. 
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The most difficult part to understand of theories with spontaneously broken gauge 
theories is, together with the Higgs potential, the Yukawa sector. Its origin is unclear, 
and probably related to the origin of the spontaneous breaking of the symmetries itself. 
It contains a large number of free parameters that must be adjusted by hand to obtain a 
realistic spectrum of particles and mixings. From the 17 parameters needed to describe the 
electro-weak theory, 13 come from the Yukawa sector. In addition, the Yukawa couplings 
contain a lot of spurious parameters that cannot be observed because they can be removed 
from the Lagrangian with an appropriate choice of phases and mixings for the fermion 
fields. Thus, a simple question such as how many parameters are needed to describe 
a theory with Yukawa couplings cannot be answered at first sight. We cannot even 
say whether the couplings are real or not and so if the theory conserves CP or not. 
This represents a great complication when one tries to understand mass matrices in the 
Standard Model. In the last years this subject has received a lot of attention (for some 
popular mass-matrix models, see for example |l|, @, [|, £§). Usually people make ansatze 
about a particular form of mass matrices and try to extract from them a reasonable 
spectrum of masses and mixings. However, completely different forms of mass matrices 
may be equivalent in their physical content. This has led to some confusion in the past. 
Moreover, some symmetries actually present in the Yukawa sector may be hidden when 
one does not work in the appropriate basis. For those reasons, an analysis based on 
symmetry should be superior and could clarify some of the issues present in particular 
forms of mass matrices. It would also allow a clear understanding of the number of 
parameters needed to describe the model. 

The standard method to calculate the number of parameters needed to describe the 
Yukawa sector of a spontaneously broken gauge theory is just by construction. After 
symmetry breaking one diagonalizes the mass matrices of all fields and after that one uses 
the freedom in the phase definition of the fields to reabsorb as many phases as possible. 
The remaining parameters in the Lagrangian are the physical parameters. As we know, 
the method works very well for the Standard Model. However, as the complication of the 
Yukawa sector increases, the method becomes intricate and can lead very easily to error; 
the question of CP-violation, for example, can then be difficult to answer. In addition, in 
many cases, phases and mixings can be moved from the gauge couplings to the Yukawa 
sector or even included in a redefinition of the fields || [|, which leads to further confusion. 

Here we propose a method, based on symmetry, to perform the counting before spon- 
taneous symmetry breaking and without performing the explicit diagonalization of the 
mass matrices. The method clarifies the structure of the manifold of the parameters that 
are needed to describe the Yukawa sector of a gauge theory and suggests a parametriza- 
tion that is useful when one wants to study the behaviour of the Yukawa couplings under 
the renormalization group. Basically the method consists in the study of the symme- 
tries of both the full and the Yukawa Lagrangians. The number of parameters needed to 
describe the Yukawa sector comes from the balance between the number of parameters 
contained in the Yukawa matrices and the number of symmetries broken by the Yukawa 
couplings. As an example, we will start with the Standard Model with the most general 
set of Yukawa couplings. We will then impose to the Standard Model a set of chiral 
symmetries that could reduce the number of parameters and see how the method can be 
used in those cases. Finally we will apply the method to more involved models based on 
the same gauge group but different particle content. 

The number of parameters needed to describe the Yukawa sector of a theory is ob- 



viously the same before and after spontaneous symmetry breaking^. Thus, we will start 
with the hadronic part of the electro-weak Lagrangian before symmetry breaking 

C = iQ L ft Ql + iuR ft u R + id R ft d R + {Q L Y u u R ip + Q L Y d d R (p + h.c.) . (1) 

Here Ql, u r and d R are the standard quark fields, left-handed doublet, u-singlet and 
d-singlet respectively. If we assume n generations they are n-column vectors in genera- 
tion space; Y u and Yd are the Yukawa couplings represented bynxn complex arbitrary 
matrices; ft = 7 M -D M where is the Standard Model electro-weak covariant derivative. 
Finally ip is the Higgs doublet and dp = iT 2 ip*. 

If Y u = Yd = 0, the Lagrangian is obviously invariant under the following chiral 
symmetries: 

Ql -> VqQ l u r -> V u u R d R -> V d d R , (2) 

where, Vq, V u , and are n x n unitary matrices acting on flavour space. The Yukawa 
couplings break explicitly these symmetries, but in a very particular way. In fact, if we 
let the Yukawa couplings transform as follows 

Y u - V Q Y u Vl Y d - V Q Y d V} , (3) 

the Lagrangian of eq. ([!]) is still invariant under the combined action of the transformations 
in eqs. (§) and (Q). It is also easy to see that not all the symmetries in eq. @ are broken by 
the Yukawa couplings. Indeed if we choose Vq = V u = V d = e ia , the Yukawa Lagrangian 
remains invariant. This is nothing else than baryon number conservation. 
Equation (^) defines an equivalence relation 

{Y U) Y d ) & (K, Y£ = (V Q Y u Vl V Q Y d Vj) . (4) 

The Lagrangians with couplings (Y u , Y d ) and couplings (Y£, Yj) are completely equivalent. 
Thus, counting how many parameters are needed to describe masses, mixings, and Yukawa 
couplings, is the same as counting how many equivalent classes there are with respect to 
the equivalence relation of eq. (f|). 

This problem has some similarity with the problem of spontaneous symmetry breaking 
by the Higgs mechanism. There, one has a group G and a representation if of Higgses 
whose vacuum expectation value breaks the group G down into the group G' C G and 
one wants to know the number of physical Higgses. As we know, the number of physical 
Higgs degrees of freedom N v . s is 

■^iPphys iV Goldstone ■ (5) 

Here N v is the number of degrees of freedom in the Higgs representation ip and Nooidstone 
is the number of Goldstone bosons that appear after spontaneous symmetry breaking. It 
is equal to the number of broken generators of G, which is the number of generators of the 
full group G minus the number of generators of the unbroken subgroup G' (N Go i dstone = 
Nq — Nqi). We have a similar situation here. A chiral symmetry G is broken explicitly 
by the Yukawa sector Y to a group G'. Then, only the broken part of G can be used to 
absorb parameters from Y. Following eq. @ we could write 

N^ ys =N Y -N G + N G , . (6) 

^or the moment, we leave aside the case of spontaneous breakdown of CP, in which some parameters 
can be moved from the Higgs potential to the Yukawa sector. 



In the Standard Model, as follows from eq. (0), the group is G = U(u)q ® U{n) u ® 
U(n)d and the Yukawa couplings Y = (Y u , Yd) are two general n x n complex matrices 
transforming under the group G as eq. (0j). Here the factors U(n) denote the full unitary 
group U(n) = SU(n) <g> U(l). As discussed previously, for the most general Yukawa 
couplings, the only subgroup of G that leaves the Yukawa couplings invariant is a U(1) R 
such that V Q = V u = V d = e ia . 

After eq. @, counting of parameters is trivial. Taking into account that an arbitrary 
complex matrix of dimension n contains n 2 moduli and n 2 phases and that a U(n) matrix 
contains n(n — l)/2 moduli and n(n + l)/2 phases, we find that Y p h ys can be expressed 
in terms of 



couplings and symmetries 


moduli 


phases 


(Y u , Y d ) 


2n 2 


2n 2 


U{n) Q ®U{n) u ®U{n) d 


-3n(n - l)/2 


-3n(n + l)/2 


U(1)b 





1 


Yphys 


2n + n(n - l)/2 


(n-2)(n-l)/2 



In particular if n = 3 we have 6 + 3 moduli and 1 phase. After symmetry breaking, some 
of the moduli are related to masses and other to mixings or remaining Yukawa couplings. 
To distinguish among them we can repeat the argumentation but using the completely 
broken theory with all the heavy and Higgs bosons decoupled. In the case of the Standard 
Model, what remains is a model with the standard fermion content but only QED (and 
QCD if included in the analysis)^. The resulting Lagrangian is 

C = iul Jp u L + HTr~ p u r + id L p d L + id R p d R + (uZM u u R + d L M d d R + h.c) . (8) 

Here p is the standard QED covariant derivative and M u and Md are the quark mass 
matrices related to the Yukawa couplings by the vacuum expectation value of the Higgs 
scalar. 

Except for the mass terms, the Lagrangian of eq. (§) is invariant under the following 
symmetries 

u L -> V UL u L u R -> V UR u R d L -> V dL d L d R -> V dR d R . (9) 

Using similar arguments as above, but now taking into account that the full Lagrangian 
has separate flavour conservation for w-type and ci-type quarks, i.e. G = U{n) uR <g) 
U{n) uL ®U{n) dR ®U{n)dL and G' = (U(l) u ) n <g> (U(l) d ) n , we find that the Lagrangian © 
gives rise to only 2n masses. Of course this is obvious since the symmetries in eq. @ are 
enough to diagonalize the two mass matrices completely. Thus, from the 2n + n(n — l)/2 
moduli present in the Yukawa sector of the Standard Model, 2n correspond to masses and 
n(n — l)/2 to mixings, as expected (see also 

Yukawa couplings related by eq. @ are in the same class of equivalence. To find 
a parametrization of the physical Yukawa sector ve have to characterize the equivalent 
classes, and this can be done by taking one element of each equivalence class. We can use 
a Vq and a V u rotation to diagonalize the matrix Y u , VqY u V^ = D u . After that, since any 
arbitrary complex matrix can be written as a Hermitian matrix times a unitary matrix, 
we can use a Vd transformation to write Yd as a positive-definite Hermitian matrix. Thus 
(Y u , Yd) — > (D u , Hd), with D u a positive-definite diagonal matrix and H u a positive-definite 
Hermitian matrix. But, D u is still invariant under transformations like D u — > KD U K\ 

2 We have to keep all the exact symmetries. Fields with different charges never mix. 
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with K a diagonal matrix of phases, while Hd is not. This means that we can use K to 
absorb phases from H u . In fact we can chose K such that the next-to-diagonal diagonal 
elements of Hd be real and positive. The only remaining symmetry is just U(l) b of baryon 
number, as expected. Thus, without loss of generality we can represent the physical 
Yukawa couplings of the Standard Model by a positive-definite diagonal matrix D u for 
the u-type quarks and a positive-definite Hermitian matrix Hd with the next-to-diagonal 
diagonal elements of Hd real and positive: 



(Y u ,Y d 



Y, 



phys 



(At, Hd 



(10) 



One can easily check that the right-hand side of eq. fll0|) contains also 2n + n(n — l)/2 
moduli and {n — l)(n — 2)/2 phases. Of course, one can further write Hd in terms of a 
Kobayashi-Maskawa-type matrix and a diagonal matrix of masses. 

Since the full Lagrangian is invariant under the combined action of transformations 
(0) and (|3p it is easy to show that the renormalization group equations for the Yukawa 
couplings must be covariant with respect to the set of transformations (§). This guarantees 
that the renormalization group equations can be written only in terms of the parameters 
in the right-hand side of eq. ([T~0|) . This way we can obtain a set of renormalization group 
equations with all the unphysical parameters removed^. 

The method we just explained for counting the physical parameters works very well 
in the Standard Model with the most general Yukawa couplings. However, in the form 
we have presented it, it is not suitable when some of the masses are zero or there are 
some symmetries among them. We want to generalize the method to those cases. We will 
impose that the only allowed reductions of parameters are those that are protected by 
some symmetry. For instance, suppose we want to make the first n — 1 ci-type quarks (in 
n generations) exactly massless. We could impose an additional chiral symmetry such as 



Ql^Ql 



Ur -> Ur 




(11) 



where U is a (n — 1) x (n — 1) general unitary matrix acting only on the first n — 1 
generations of cLr quarks. Clearly the only way to make the Yukawa Lagrangian invariant 
under this symmetry is for 



/ 




Y u = arbitrary Y d 



Qtl \ 

a 2 







a, 



(12) 



/ 



Then, the Yukawa Lagrangian is only invariant under U(n — 1) <£> U(1)b- Naive counting 
starting from the full group G = U(n)g g) U{n) u g) U{n) d does not work in this case, 
because the couplings in eq. (|T2|) are not a representation of the full group G. Said in 
another way, general G transformations do no preserve the form of the Yukawa matrices 
in eq. (0) and, thus, we cannot use the full group to absorb parameters from the Yukawa 
sector. It is not difficult to see that the only subgroup of G that preserves the form of the 
Yukawa couplings in eq. ( |I2"D is G = U (ti)q <g) U {n) u £g> U{n — \)®U{1)b- Now we can use 
G to absorb parameters from the Yukawa sector and the counting of physical parameters 



3 A similar approach has been followed by [|[ |To[ [TT| [l2| 
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comes as follows 



couplings and symmetries 


moduli 


phases 


(Y u ,Y d ) 


n 2 + n 


n 2 + n 


U{n) Q ®U{n) u ®U{n- l) d ®U{l) B 


-(n-l)(3n-2)/2 


-n{3n+ l)/2 - 1 


U{n-l)®U{l) B 


(n - l)(n - l)/2 


(n - l)n/2 + 1 


Yphys 


n + 1 + (n - 1) 






• (13) 



The model leads to n massive w-type quarks, 1 massive d-type quark and n — 1 mixings. 
There is no CP-violating phase. This result can be explicitly checked by full diagonaliza- 
tion of the Yukawa sector. 

Another extreme case can be obtained by imposing an axial U(l) symmetry on the 
d- quark: 



Ql — > Ql u r —> u r d iR — > d iR i = 2, 
Then, the Yukawa couplings take the following form 



n 



d 



1R 



e ia d 



1R 



(14) 



Y u = arbitrary 



Y d 



( a 1;1 
a 2 ,i 



\ a n ,i 



«l,n-l \ 



a, 



(15) 



It is easy to see that the group that preserves the form of these matrices is G = U(n)q <g> 
U (n) u <S> U'(n — l)d ®U(1)b with the U'(n — 1) now acting on the 2, 3, • • • , n generations 
of d R quarks, while the only symmetries of the full Lagrangian are 11(1)^ ®C/(1)b- Again 
we can do the counting 



couplings and symmetries 


moduli 


phases 


(Y u ,Y d ) 


n 2 + n(n — 1) 


n 2 + n(n — 1) 


U{n) Q ®U(n) u ® U'{n - l) d ®U(l) B 


-(n-l)(3n-2)/2 


-n(3n + l)/2 - 1 


U{l) dl ®U{l) B 
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y 

1 phys 


n + (n-l) + n(n - l)/2 


(n - l)(n - 2)/2 



(16) 

This model just gives a massless <i-quark. Phases and mixings are exactly the same as 
in the Standard Model for any number of generations. This result is interesting because 
the so-called m u = solution of the Strong CP Problem (see for example [fT3|, [14], 15| and 
references therein) is based in a symmetry like ours. In our discussion we have obviated 
the fact that chiral symmetries have in general anomalies. In the Standard Model, be- 
cause of QCD instanton effects and because the chiral symmetries are anomalous, there 
is an additional phase that cannot be removed. Chiral transformations move it from the 
Yukawa sector to the QCD sector. This represents a problem because it gives a too large 
contribution to the electric dipole moment of the neutron unless it is made unnaturally 
small. If for instance m u = 0, there is an additional chiral symmetry (similar to the one 
we have just studied) and the extra phase can be rotated away. However, imposing an 
additional chiral symmetry is a very strong constraint on the mass matrices [see eq. (fI5|)]. 
One has to redo all parameter counting and be sure that, there still is in the model an 
observable CP-violating phase that could account for all observed CP-violating effects. 
This is just what we did in our example. 



The examples discussed above show that the use of symmetry is a natural and basis- 
independent way of choosing models for mass matrices. Of course, most of the symmetries 
lead to some massless quarks and/or some vanishing mixings, but the symmetries could 
be softly broken in order to generate a realistic pattern of masses and mixings. We 
find it interesting to explore and classify all possible symmetries that can be imposed 
on the Yukawa sector of the Standard Model according to the spectrum of masses and 
mixings they generate. This is nothing but the classification of all subgroups of U {ti)q ® 
U(n) u ® U{n)d according to the spectrum they produce. We think that this classification 
is essential for any mass-matrix modelling. 

Let us apply our method of parameter counting to a more complicated example. Let 
us take an extension of the Standard Model including right-handed neutrinos with a 
Majorana mass term. The leptonic part of the Lagrangian is 

C = ih ft h + ivr ft vr + ieii ft e R + (hY u v R <p + hY e e R fi + ^%Mv R + h.c)) (17) 

where £l is the standard leptonic doublet, e R is the right-handed charged lepton singlet, 
v R is the right-handed neutrino, and v c R = Cv^ . The Yukawa couplings are the usual Y e 
that give masses to the charged leptons, Y v that give a Dirac mass term to neutrinos and 
M is a complex symmetric Majorana mass matrix for right-handed neutrinos. It can be 
put in by hand, as we did, because right-handed neutrinos are singlet, however, it could 
be generated through the vacuum expectation value of some extra scalar singlet. The 



model just described is nothing else than the see-saw model [16|, [T7|, 18, |1^| of neutrino 



masses. IfY u = Y e = M = the Lagrangian is invariant under the following symmetries 

h^Vdi vr^V v v r e R ^V e e R , (18) 

Ve, V v , and V e are unitary matrices of dimension n for n generations of leptons. As in the 
Standard Model, the Yukawa couplings break these symmetries explicitly. However, they 
are broken in such a way that if we let the Yukawa couplings and the Majorana mass 
transform in the following way 

n-W£ Y u ^V e Y u Vj M^V u *MVj, (19) 

the Lagrangian remains invariant. At difference with the Standard Model, we see that 
the Yukawa couplings and the Majorana mass term completely break all the invar iances 
of eq. (|T8"|), including lepton- number conservation. 

Equation (^) defines the following equivalence relation 

(Y e , Y v , M) (Y e ', Yl, M') = (V t Y e Vj, V e Y v V*, V*MVj) . (20) 

Again the number of physical parameters needed to describe the model is given by the 
difference between the number of parameters contained in the Yukawa coupling and the 
number of generators broken by the couplings. 

Now the group of symmetries of the kinetic and gauge part is G = U(n)e <g) U(n) e <g> 
U{n) v and there is no unbroken subgroup. Counting of parameters is easy. Taking into 
account that an arbitrary complex matrix of dimension n contains n 2 moduli and n 2 
phases, that a symmetric complex matrix contains n(n + l)/2 moduli and n(n + l)/2 
phases, and that a U{n) matrix contains n{n — l)/2 moduli and n{n + l)/2 phases, we 
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find that Y p h ys can be written in terms of 



couplings and symmetries 


moduli 


phases 


(Y e ,Y„,M) 


2n 2 + n(n + l)/2 


n 2 + n(n + l)/2 


U(n) e ®U(n) e ®U(n) v 


-3n(n - l)/2 


-3n(n + l)/2 
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1 phys 


3n + n(n — 1) 


n(n — 1) 



We obtain that even in the two-generation case there are two observable phases that 
could violate CP. As we did in the Standard Model, to know how many of the physical 
parameters are related to masses and how many to mixings, we have to decouple all kinds 
of massive bosons. Then, all the mass matrices can be diagonalized, leading in general 
to n massive charged leptons and 2n Majorana neutrinos. Thus, 3n of the moduli are 
masses, and the rest correspond to mixings. 

To find a parametrization of the Y p h ys , we can use the equivalence relation in eq. (|20|) 
to reduce the number of parameters in the Yukawa couplings. Using the fact that an 
arbitrary complex matrix can be written in a unique way as a unitary matrix times a 
positive-definite Hermitian matrixF] and that a complex symmetric matrix can be written 
in a unique way in terms of a positive-definite diagonal matrix and a unitary matrix 
S = U T DU we can write the Yukawa couplings in the following form 

(Y e , Y v , M) -> Y phys = (H e , H v , D m ) . (22) 

Here H e and H v are positive-definite Hermitian matrices and D m is a positive-definite 
diagonal matrix. As we have exhausted all the symmetries in the group G to reduce the 
Yukawa couplings to this form and there is no subgroup of G that leaves Y p h ys invariant, 
the set of parameters in the right-hand side of eq. ( [22] ) is the physical set of parameters 
needed to describe the theory. This means that we could start from the beginning without 
loss of generality by choosing Y v and Y e as positive-definite Hermitian matrices and M 
as a positive definite diagonal matrix and all the parameters in those matrices would be 
observable as can be easily checked by counting the number of parameters contained in 
those couplings. 

Finally to give a flavour of how these techniques can be used to study more exotic 
models, we chose the model of refs. 0, [H], |22|, |23| with the following Lagrangian, 

C = i£ L ft £ L + ivji ft vr + ieR ft e R + isZ ft s L 

+ (l2Y u v R Lp + T L Y e e R ^ + -Vr-M s s l + h.c)j , (23) 

where sl is a singlet that carries the same lepton-number as the leptons and M s is an 
arbitrary complex matrix. The rest of the notation is the same as in the previous example. 
In the Lagrangian of eq. (|23|) , total lepton-number has been imposed as a global symmetry. 
There are thus no Majorana mass terms for any of the singlet fermions. The group of 
invariances of the non- Yukawa part of the Lagrangian is U(n)t <g) U(n) e <g> U(n) v <8> U(n) s , 
where U(n) s is the new invariance related to the field sl- The only symmetry of the 
Yukawa couplings is just U(l)i ep of lepton-number conservation. Doing now the usual 

4 We will assume that the matrices have non-zero determinant and there are no degeneracies. If there 
is some degeneracy or zero eigenvalues, this must be imposed with some symmetry and would change the 
counting of parameters. 
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counting we get 



couplings and symmetries 


moduli 


phases 


(Y e ,Y u ,M s ) 


3n 2 


3n 2 


U (n) t ® U (n) e <8> C/ (7i)„ <8> C/ (n) a 


-4n(n - l)/2 


-4n(n + l)/2 


«7(l)iep 





1 




2n + n 2 


(n - l) 2 



(24) 



Once the gauge symmetry is broken, the Yukawa couplings generate mass terms for the 
leptons. It is easy to see that, as a consequence of lepton-number conservation, all neu- 
trinos must be Dirac neutrinos or just massless. Then, from the full mass matrix of 
neutrinos, n massless neutrinos and n massive Dirac neutrinos arise, hence only n masses 
come from the neutrino mass matrix. The other n masses come from the charged lepton 
sector. Thus, from all the moduli, 2n of them correspond to masses and the rest n 2 cor- 
respond to mixings. In addition there are On — l) 2 CP- violating phases that cannot be 



removed. This result is in complete agreement with the results obtained in [23] after full 
diagonalization of the mass matrices. 

We have shown that a study of the chiral symmetries of the kinetic and gauge parts 
of the Lagrangian and the Yukawa part of the Lagrangian before spontaneous symmetry 
breaking allows us to compute, in a straightforward way, how many observable parameters 
come from the Yukawa couplings in a general gauge theory. Basically, we have to know 
which subgroup of chiral symmetries of the kinetic and gauge parts of the Lagrangian 
preserves the form of the Yukawa couplings and which group leaves them invariant. The 
number of parameters that can be absorbed from the Yukawa sector is just the difference 
between the parameters needed to describe those groups. Then, the number of observable 
parameters comes from the balance between the Yukawa couplings and the symmetries 
broken by them. The method is useful when the diagonalization of the full set of mass 
matrices becomes complicated. On the other hand it clarifies the analysis of observable 
parameters. We explain the technique with several examples based on the Standard 
Model, with Yukawa couplings obeying various chiral symmetries, and also with some 
extensions of the Standard Model with enlarged fermion sector. 

The study of the mass matrices and the gauge couplings is not general enough, because 
some parameters such as phases, mixings, etc, can be moved from charged currents to 
Yukawa couplings or even be included in a field definition || |6|. Only the analysis of 
the full Lagrangian, and this is better done before symmetry breaking, can report the 
number of observable parameters in a non-ambiguous way. This point of view suggests 
a parametrization of Yukawa couplings in terms of invariants with respect to the chiral 
symmetries of the kinetic and gauge parts of the Lagrangian that can be useful in the study 
of the behaviour of Yukawa couplings under the renormalization group. We comment on 
this parametrization in the case of the Standard Model and the see-saw model for neutrino 
masses. 

We find that a complete analysis and classification of all the chiral symmetries that 
can be imposed on the Yukawa sector, according to which spectrum of masses and mixings 
they originate is still missing. We think this analysis would be a good starting point for 
mass matrix modeling. 

Finally we want to remark that we only considered symmetries acting on fermions. 
The analysis could be extended to symmetries involving fermions and scalars as well, like 
the Peccei-Quinn f25 , 26] symmetry or horizontal symmetries |27|, p8 , 
think that it can be extended to the spontaneous CP-violation case. 



2|, |3 . We also 
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